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The long-range proximity effect in superconductor/ferromagnet (S/F) hybrid nano- structures is observed if 
singlet Cooper pairs from the superconductor are converted into triplet pairs which can diffuse into the fer- 
romagnet over large distances. It is commonly believed that this happens only in the presence of magnetic 
inhomogeneities. We show that there are other sources of the long-range triplet component (LRTC) of the con- 
densate and establish general conditions for their occurrence. As a prototypical example we consider first a 
system where the exchange field and spin-orbit coupling can be treated as time and space components of an 
effective SU(2) potential. We derive a SU(2) covariant diffusive equation for the condensate and demonstrate 
that an effective SU(2) electric field is responsible for the long-range proximity effect. Finally, we extend our 
analysis to a generic ferromagnet and establish a universal condition for the LRTC. Our results open a new 
avenue in the search for such correlations in S/F structures and make a hitherto unknown connection between 
the LRTC and Yang-Mills electrostatics. 



PACS numbers: 74.45,+c, 74.78.Na, 75.70.Tj 

The odd triplet superconductivity in superconductor- 
ferromagnet (S/F) structures has been intensively studied, 
both theoretical and experimentally, since its prediction in 
2001 dEI. In that context it is of crucial interest to under- 
stand the process of converting the singlet Coppers pairs from 
the superconductor into triplet pairs of electrons with equal 
spins in the ferromagnet. Apart from its interest for funda- 
mental research, the study of triplet superconducting correla- 
tions might find useful applications for spintronics 0. 

It is well established that triplet pairs, once created, diffuse 
into the ferromagnetic materials over distances much larger 
than the singlet ones. This leads to the long-range proximity 
effect which explains the observation of Josephson currents 
through S/F/S junctions over large distances PUTOl and the 
long-range propagation in S/F structures of superconducting 
correlations ifTTIl . The usual theoretical interpretation of these 
experiments assumes that the singlet- triplet conversion is me- 
diated by a magnetic inhomogeneity in the vicinity of the S/F 
interface EH. This can be caused by a domain wall [2], a spin- 
active S/F interface |' 13], or by a multilayered ferromagnetic 
structure with different magnetic orientations lfT4ll . 

Formally, the existence of the long-range triplet component 
(LRTC) can be inferred by inspecting the spin structure of 
the quasiclassical condensate (anomalous) Green's function, 
/ = f s + fi<7. Here f s is the amplitude of the singlet com- 
ponent, and the vector f t describes the triplet correlations (d 
is the vector of Pauli matrices). The component of the vec- 
tor f t parallel to the magnetization describes the triplet state 
with zero spin-projection (i. e., | t,|) + I Lt)) fI51H51 . The 
LRTC corresponding to the pairs with spin-projections ±1 
only exists if f t is non-collinear with the magnetization direc- 



tion, which is the case for certain magnetic inhomogeneities 
[03HI21- Indeed, it is commonly believed that the only way to 
generate the LRTC in S/F hybrids is by means of creating a 
non-homogeneous magnetic configuration. 

In this Letter we study the superconducting LRTC in S/F 
structures from a more general perspective and demonstrate 
that, besides an inhomogeneous magnetization, there are other 
sources of long-range triplet correlations. In particular, the 
momentum dependence of an effective exchange field,which 
can be attributed to the spin-orbit (SO) coupling naturally gen- 
erates LRTC, provided certain conditions are fulfilled. 

We show that the physical mechanism of the singlet-triplet 
conversion can be linked to the local SU(2) invariance of mag- 
netized systems with SO interaction QjOEni- To reveal this 
physics we first analyze a prototypical example of a spin- 
dependent field consisting of a momentum independent Zee- 
man term and a linear in momentum SO coupling. These two 
contributions act as the time and space components of the 
SU(2) gauge potential, which ensures that that the Zeeman 
and SO fields can enter physical quantities only in gauge co- 
variant combinations, such as the SU(2) field strength l2T1i . 
We derive a SU(2) covariant Usadel equation and identify 
the SU(2) electric field as a key object responsible for the 
long-range triplet proximity effect and for the existence of the 
LRTC in S/F structures. By solving the Usadel equation for 
a lateral S/F junction we demonstrate that even in the case of 
a uniform exchange field, the LRTC is present in the system. 
This gives a plausible explanation for the experimental obser- 
vation of long-range Josephson effects in S/F/S junctions. 

In the second part, we generalize our results to systems 
with possibly anisotropic Fermi surfaces and generic spin- 
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dependent fields with an arbitrary momentum dependence. 
We derive the general quasiclassical equations for the anoma- 
lous Green's function, which allows us to establish universal 
conditions for the creation of triplet long-range superconduc- 
tivity in diffusive S-F hybrids. Our theory gives an unified 
physical picture of the LRTC and opens up a possibility to 
predict long-range proximity effects in S/F structures from the 
knowledge of the electronic structure at the Fermi level. 

We consider a hybrid structure consisting of a conventional 
s-wave superconductor (S) with the order parameter A in con- 
tact with a ferromagnet (F) with a SO coupling. Let us first 
assume that SO effects can be modeled by a generic linear 
in momentum form H$o = ^{Pj^j}^ where Aj = A a -G a 
are the components of a 2 x 2 matrix valued vector that 
parametrizes the SO coupling. In this case the Hamiltonian 
in the F-region can be represented in the form 



Here V k > = d^- —i[Ak,-] is the covariant gradient operator. 
At the S/F boundary with normal vector N , we use the 
Kuprianov-Lukichev boundary conditions | 25 ] which take the 
form 



tf*V*/| 7 = -y/ A , 



(5) 



-A + Vi 



impi 



(1) 



where /a = A/V&) 2 + A 2 is the anomalous Green's function 
in the S -region. Hence in the covariant formalism the usual 
gradients are replaces by the covariant ones, and this is the 
only place where the SO coupling enters the theory. Eqs. Q, 
§5§ are manifestly gauge covariant and their structure is very 
appealing physically. In fact, they can be written immediately 
by using only the gauge symmetry arguments [27] . 

Next, we write / as the sum of the singlet f s and triplet 
ft contributions, splitting out of f the part parallel to the ex- 
change field Aq: 



where Ao — A^a a = h a a a is the exchange field, and Vi mp is 
the impurity scattering term. The SO coupling and the Zee- 
man term enter the problem as the space and time components 
of the SU(2) gauge potential, which implies the SU(2) gauge 
invariance. The Hamiltonian remains unchanged under any 
local SU(2) rotation with a matrix U (r) supplemented with 
and the gauge transformation of the potentials, 



f = fs+ft =fs+A f?+f t L . 



(6) 



For any matrix M we have defined M 1 - = 
\[A ,[Ao,M]]/\A \ 2 , where \A \ = \jA a A a is the am- 
plitude of the exchange field. 

Trace of Eqs. ([4]), ([5} gives the equations for the singlet 
amplitude f s coupled to the parallel triplet amplitude f} : 

.sgnco 



Aj^UAjU- l -i(djU)U-\ Ao^UAoU- 1 . (2) ^ fs ~ Kfs - 2i—^-\A \ z f t 11 = 0, NAf^ = -yf A (7) 



Since we are interested in equilibrium quantities we work 
with the Matsubara 4x4 matrix (in the Nambuxspin space) 
Green's function 80(1*1,1*2) at the discrete frequencies (0 = 
7iT(2n +1). To keep track of exact SU(2) gauge symmetry 
we employ a technique developed in the context of quark- 
gluon kinetics [ 22 ] and used recently to describe spin dynam- 
ics in semiconductors (23). Namely, we introduce the covari- 
ant Green's function as follows 



where V 2 is the usual Laplace operator, and k* 2 = 2\co\/D. 
The traceless part of Eqs. ([?]), ([5]), can be rearranged as follows 



a 



D 



Sfl,(ri,r2)=W(R,ri)g fl ,(n,r 2 )W(r2,R) 



(3) 



where W(R,r\) and W(r2,R) are the Wilson link oper- 
ators which "covariantly connect" the arguments of the 
Green's function to the "center-of-mass" coordinate R = 
Yl \ Yl 1 26 ]. The advantage of the Green's function ([3J) is that 
its Wigner transform, and thus the corresponding quasiclas- 
sical Green's function |(n,R), transform locally covariantly 
under a nonuniform SU(2) rotation, i. e. g 1— >• UgU~ l . By 
using the method of Ref. l22l we can derive the equation of 
motion for 5co and then proceed further to the quasiclassical 
limit and eventually to the diffusive Usadel equation [24]. The 
detailed derivation of the full quasiclassical equation will be 
presented elsewhere. Here we only show the final linearized 
Usadel equation in the ferromagnet for the covariant anoma- 
lous function / 



N k (A d k fj l + Vkf 1 + $ w f} ) | 7 = 0. 
Here is the SU(2) electric field defined as 

$k0 = dkAo - i[A k ,Ao] . 



0, (8) 
(9) 

(10) 



Inspection of Eqs. (8|9 ) shows that 3^ is the key object for 
the long range proximity effect: If 3^ = the homogeneity 
of the boundary condition ^ implies that f^ =0. In other 
words, no transverse triplet is generated in the absence of the 
SU(2) electric field. 3^ vanishes in the case of a uniform 
magnetization and if the SO and Zeeman terms commute. 

If 3>o is small Eqs. ^ and ^ can be treated perturbatively. 
The leading contribution is given by the terms proportional to 
Ao, i. e., by the first line in Eq. ([8]) and by the first term in 

Eq. d9l). This yields the well known equations for the short- 

II 

range triplet component f t u coupled to the singlet one f s , 



DV k (V k f)-2\co\f-isgnco{A J} = . 



(4) 



v 2 /; 1 - Kifi 



-2^f s 
D 



0, N k d k ff\. = 0. (11) 
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Equations for the lowest in 3^ correction to f t come from the 
second line in Eq. and the last two terms in Eq. d9l): 



v 2 // 



(12) 
(13) 



In the first terms in the left hand sides we replaced with 
as the difference of these operators gives higher order correc- 
tions compared to those determined by the terms ~ f} . 

Equations ([7]), ([TT])-([T3]) provide a complete description of 
S/F structures, which clearly demonstrates a common physical 
role of the SO coupling and inhomogeneous magnetization in 
the problem of the singlet- triplet conversion. As expected on 
general grounds, they appear in the theory in form of a single 
gauge covariant object - the SU(2) electric field 3>o Eq. (10) 
entering the "source part" of Eq. ( 12 ) for f t L . 

It is easy to see that the well known generation of LRTC 
by magnetic inhomogeneities follows naturally from our co- 
variant formulation. Consider for example a transversal mul- 
tilayer S/F structure shown in Fig la in the absence of SO 
coupling. In this case a nonzero 3^ is solely due to inhomo- 
geneity of the exchange field Ap. Assume that Aq has only an 
in-plane component (eventually rotating). Then the first term 
in the r.h.s. of Eq. ( [TO] ) generates the LRTC for a Bloch do- 
main wall parallel to the interface [2], while the second term 
~ Vfc3>o is responsible for the LRTC in the presence of a fi- 
nite Neel wall perpendicular to the interface plane [30]. It is 
interesting to note that the covariant derivative Vk&kO of the 
non-Abelian electric field is exactly the right hand side of the 
Yang-Mills electrostatic equation. The general gauge symme- 
try arguments of Ref . I2H (used there to uncover the nature of 
the equilibrium spin currents) show that V^3^ is proportional 
of the magnetization induced in the F-region by a non-uniform 
exchange field and/or SO coupling. This reveals the nature of 



the second term in Eq. ( [12] ) and provides an interesting con- 
nection between the generation of LRTC at the edges of Neel 
domain walls, and the Yang-Mills electrostatics. 

We now analyze a SO-generated LRTC in the case of 
an uniform magnetization. A special type of SO coupling 
in Eq. ([T]) should naturally occur in the vicinity of hetero- 
interfaces where inversion asymmetry exists OH . Hence we 
concentrate on this particular situation when the SU(2) vector 
potential is localized around the S/F interface and has only in- 
plane components A x and A y . This implies that only 3*0 and 
3^0 are nonzero. 

As a first example we consider a transversal S/F structure 
(Fig. la) assuming for definiteness a Rashba-Dresselhaus SO 



term with A x = j36 x 



aa y and A y 



(ccd x — pG y ), where a 



and j3 are the Rashba and Dresselhaus constants. We assume 
a contestant in-plane magnetization. Thus, only the second 
term in the r.h.s of Eq. ( p"2| ) serves as source for the LRTC. 
One can easily show that this term is non-zero only if a/3 ^ 
and Aq ^ A y Q . This in particular means that a pure Rashba 
or Dresselhaus SO coupling does not induce the LRTC in a 
transversal geometry with an in-plane magnetization. 



Lateral S/F structures are more favorable for the existence 
of LRTC. Consider the structure shown in Fig. lc. Assuming 
for simplicity, but without loss of generality, that the F film is 
thin enough, we integrate Eqs. ([7]), ([TT])-([T3]) over z-direction, 
and obtain the following set of ID equations 



dx fs — K a)fs - 

ft ~ K Q)ft^ 

d 2 Jt L -*lft 



-2i 



sgnco - 2 



D 

2/sgnco 

D 



\A \ z fl l = -0(-x)ff A , (14) 



/* = 0, 



-26(-x)? x0 d x f}\ 



(15) 
(16) 



where y and 3^ are effective values averaged over the thick- 
ness. The boundary conditions atx = are the continuity of all 
functions, and the continuity of d x f s , d x f} , and d x f^~ + 3x0 ft • 
This boundary problem can be solved straightforwardly. In 
the interesting limiting case of the exchange field |.Ao| much 
larger than the superconducting gap A the condensate function 
at x > (cf. Fig. lc) is given by the expression 



/(*)= L C( * e 
a=l,2 



-^Ce- K °> x , (17) 



where C h2 * -y/ A £ 2 Kf,i/(4|A)| 2 ), K?, 2 = ±2isgnfi)|yio|/A 
with Re/Ci 5 2 > and C « -(3//2)7/ A Z)sgna)/(|yL | 2 K- a) ). All 
components decay away from the edge plane x = 0. The 
singlet and parallel to Aq triplet components [the first term 
in Eq. (17)] decay over the short magnetic distance ^/D/2h, 
while the triplet component perpendicular to Aq [the second 
term in Eq. ( [17] )] decays over a larger length of the order of 
y/D/2T confirming its long-range character. Figure lc shows 
tha the LRTC decays in both directions from the inhomogene- 
ity at x = 0, which looks very similar to the LRTC generated 
at the edge of a Neel domain wall (30). This similarity is not 
accidental. In fact, in the particular case of A y = our system 
is gauge-equivalent to the Neel wall with an edge at x = 0. 

We now consider a lateral structure with two S -electrodes 
separated by a distance L (see Fig. lb). If L ^> ^D/2h the 
Josephson coupling is only mediated by the LRTC, and the 
critical current is given by 



Ir = 



tr(^ ) 2 T^ K coC 2 (co n )e- K - L , (18) 



(O n 



where S and Op are the cross section and conductivity of the 
F-region. From this equation one concludes that a finite SU(2) 
electric field with a component perpendicular to the magneti- 
zation is the source of the long-range Josephson effect. No- 
tice that the lateral geometry shown in Fig. lb is equivalent 
to the one explored in the experiment of Ref. (4). One argues 
in that case that the long-range Josephson current is due to a 
SO coupling at the S/F interfaces. Our theory can also give an 
plausible explanation for the long-range effect observed in Co 
nanowires (8), assuming a Rashba-type SO coupling due to 
the quasi- ID geometry, which is non-collinear with the mag- 
netization of the wire. 
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FIG. 1 . Different geometers considered in the text, (a) The transver- 
sal S/F/S junction, (b) The lateral S-F structure and space depen- 
dence of F = Y.(D \f\ 2 f° r me singlet and LRTC normalized with the 
asymptotic value at x = — oo. Here £o = y/D/A, h = 10A, r = 0.05 A 
and A is the order parameter in the superconductor, (c) The lateral 
Josephson junction 



We finally generalize our results to ferromagnets with a 
generic momentum dependent effective exchange field. Our 
starting point is the following Hamiltonian in the F-region 



H 



(19) 



Here § (p) is the spin-independent part of the quasiparticle en- 
ergy. The spin-dependent contribution is written as the sum of 
an even, h{— p) = h(p), and an odd, b(— p) = — £(p), in mo- 
mentum parts. The even part h(p) describes the Zeeman-type 
exchange term, which is, in general, momentum dependent in 
realistic systems. The odd part S(p) corresponds to a generic 
SO interaction which preserves the time reversal symmetry. 

Following the standard procedure (see for example 
Ref. l28l ), one derives the Eilenberger equation [29] for a 
generic spin-dependent Hamiltonian 



v£(n)3tf+[G>T 3 ,i 



-i[h(n)z 3 U(n),g] = -^[(g),g] (20) 



where n is a unit vector pointing in the direction of momen- 
tum, (n) are components of the Fermi velocity, t is the im- 
purity scattering time, and (...) denotes averaging over n. This 
equation allows for a general anisotropy with different veloci- 
ties and spin splittings at different points on the Fermi surface. 
We focus here on the diffusive limit, in which T -1 determines 



the largest energy scale in Eq. (20). In this case Eq. (20) re- 
duces to the Usadel equation for the angle- averaged Green's 
function l24ll . For a general anisotropic ferromagnet the lin- 
earized Usadel equation takes the form 

D kj d k djf-2cof-isgnco{(h) J}-2iT[{v F k b)Af] ~ (21) 
-ix [d k (vfb)J]-T([b, [bJ]])-T({Sh,{8hJ}})=0. 



where = T(v£vf ) is the tensor of diffusion coefficients, 
and Sh(n) = h(n) — (h) is the variation of the Zeeman field at 
the Fermi surface. 



With the exception of the last term containing Sh, Eqs. (21 
and the gauge covariant Usadel equation, Eq J4]), are struc 
turally equivalent. The first three terms in Eq. ( [2l] ) correspond 
to Eq. ([4]) without the commutators coming from the covariant 
gradients. These terms in Eq. ([?]) correspond to the commuta- 
tors in Eq. |2l| involving the SO field b. Because of this simi- 
larity our analysis of Eq. ([4]) is directly applicable to Eq. ([21 



Hence we conclude that the last four terms in Eq. $2A) serve as 
a the source for the LRTC. More precisely, the LRTC is gener- 
ated if any of these terms has a finite component perpendicular 
to the exchange field (h) averaged over the Fermi surface. We 
can draw a remarkable conclusion from this result: From the 
knowledge of the electronic properties at the Fermi level of 
S/F systems, namely, from §(p), ^(p), and S(p) in Eq. (19) 
one can easily infer whether or not the LRTC would exist in 
the hybrid structure. Moreover, Eq. plj ) is quite general and 
can be used in a broad context of problems involving super- 
conductivity and spin-fields. 

In conclusion, we have presented a general description 
of the long-range triplet superconductivity in S/F structures. 
Starting from the linear in momentum SO coupling we have 
developed the SU(2) covariant theory describing the diffusion 
of the condensate and identified the SU(2) electric field as a 
physical source for the LRTC. We have also considered the 
case of an arbitrary momentum dependence of the spin-fields 
and derived a useful equation from which, by knowledge of 
the electronic structure of the ferromagnet and the interfaces, 
one can directly predict whether the LRTC is generated or not. 
Our results not only unify in an elegant way all models de- 
scribing the long-range proximity effect in S/F structures, but 
also predict new sources for the singlet-triplet conversion and 
provide a useful tool in the search for triplet superconducting 
correlations. 
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